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It is known that the presence of background magnetic field in cosmic plasma distorts the acoustic 
peaks in CMBR. This primarily results from different types of waves in the plasma with velocities 
depending on the angle between the magnetic field and the wave vector. We consider the conse- 
quences of these effects in relativistic heavy-ion collisions where very strong magnetic fields arise 
during early stages of the plasma evolution. We show that fiow coefficients can be significantly 
affected by these effects when the magnetic field remains strong during early stages due to strong 
' induced fields in the conducting plasma. In particular, the presence of magnetic field can lead to 

enhancement in the elliptic fiow coefficient V2- 
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PACS numbers: 52.35.Bj, 25.75.Ld, 12.38.Mh 

It has been recently shown by us that a deep correspondence exists between the physics of inflationary density 
fluctuations in the early universe which result in the CMBR acoustic peaks and the physics of flow in relativistic 
heavy- ion collision experiments (RHICE) f]]. We further showed that important features such as acoustic peaks and 
_ suppression of super horizon fluctuations may be present in a plot of root-mean square values of flow coefficients w™*. 
II ] The possibility of presence of this suppression of superhorizon fluctuations has been recently pointed out by Sorensen 
Qh' in R.HIC data [2|. There have been several works now which discuss flow fluctuations for a large range of values of n 
^ ' [sj . This non-trivial connection between the superhorizon fluctuations of inflationary universe and similar fluctuations 
in RHICE is now discussed in several works and its various consequences are explored (see, e.g. [4|). We mention 
that such a connection between physics of RHICE and that of inflationary universe was never anticipated earlier, 
and indeed, at first sight, it looks surprising that a concept like superhorizon fluctuation which arises from highly 
^ ' non-trivial, superluminal expansion phase of the very early universe could have any relevance for relativistic heavy- 
ion collision experiments in laboratory. Such superhorizon fluctuations in RHICE originate from the fact that in the 
center of mass frame the thermalization (and any local homogenization) happens rather quickly, within about 1 fm. 
' Initial parton energy density distribution from HIJING shows that transverse fluctuations (arising from localization 
• . of partons inside initial nucleons, and from the fluctuations in nucleon coordinates) with wavelengths significantly 

' larger than 1 fm are necessarily present at the time 1 fm even in central collisions. 
I , It was also emphasized in [Tl that various analysis tools of CMBR anisotropics can be effectively utilized for RHICE. 
' It was proposed in [IJ that instead of focusing on the average values of the flow coefficients Vn for small values of 
L" ! Hi! one should calculate root-mean square values of the flow coefficients v™" for a large range of n upto 30-40. 
. 5^ ■ Further, these calculations should be performed in a lab fixed frame, which eliminates the difficulties associated with 
determination of event plane for conventional elliptic flow analysis for non-central collisions. It was shown in [p that 
$_i : a plot of values of v™^ vs. n can be used for directly probing various flow coefficients, in particular, the elliptic ffow 
_C^j fQj- non-central collisions [6|, without any need for the determination of event-plane. 

In this paper we continue to explore this fertile connection between CMBR physics and RHICE. It was shown in 
ref. [7] that the presence of background magnetic field in cosmic plasma can distort the acoustic peaks in CMBR. 
This happens due to the presence of different waves in the plasma, with velocities depending on the angle between 
the magnetic field and the wave vector. Presence of very strong magnetic flelds in the plasma (of order 10^^ Tesla) 
during early stages in RHICE has been explored extensively recently in connection with the exciting possibility of 
observing CP violation effects [8 -10]. An important effect of the presence of such strong magnetic fields in the plasma 
will be to lead to strong variations in velocities of different types of waves in the plasma. In particular the velocity 
varies with the angle between the wave vector and the direction of the magnetic field. It is thus obvious that this 
may qualitatively affect the development of anisotropic flow. We argue that the flow coefficients can be significantly 
affected by these effects. In particular, the presence of magnetic field can lead to enhancement in the elliptic ffow 
coefficient V2 by almost 30 %. (Note, we use V2 with the present definitions to denote the elliptic flow even though we 
do not adopt the conventional usage of the eccentricity for defining the corresponding spatial anisotropy.) This raises 
the interesting possibility whether a larger value of rj/s can be accommodated by RHIC data when these effects are 
incorporated using full magnetohydrodynamical simulations. 

An important issue here is the time scale over which the magnetic fleld remains strong. The magnetic fleld arising 
from the valence charges of the initial nuclei peaks to strong values for a very short time, essentially the passing time 
of the Lorentz contracted nuclei (^ 0.2 fm for RHIC energies). Subsequently it rapidly decays (^ t~^) [3]. In such a 
situation the effect of magnetic fleld on flow coefficient will be suppressed as time scale for the development of flow 
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is several fm. Even at LHC, where expected magnetic field is more than an order of magnitude larger than that at 
RHIC, the effect of this initial pulse of magnetic field on flow coefficients may not be very significant, though due to 
various uncertainties, a full magnetohydrodynamics simulation is needed to investigate these issues. 

However, it has been recently pointed out fll'l that magnetic fields of similar magnitude, as the peak value of the 
external field, can arise from induced currents due to rapidly decreasing external field. Further, in the quasi-static 
approximation, the magnetic field satisfies a diffusion equation with the diffusion constant equal to l/((T/i) where /i 
is the magnetic permeability and a is the electrical conductivity [l^. It is argued in ref. [ll| that due to this the 
magnetic field may survive for much longer time, and can lead to interesting effects. 

We will take ~ 1 (refs.[lll However, the value of electrical conductivity we take is smaller than what is 

taken in ref. fll'l. We use a ~ O.ST (= 0.3 fm"! for T ~ 200 MeV) from refs.fll]. The time scale t over which the 
magnetic field remains essentially constant [llj over length scale L is, r ~ For L = 5 — 6 fm, we get r ~ 2 — 3 



fm. For higher temperatures a will be larger increasing the value of t. (See, also, ref. [ij] for effects of leptons on ct, 
though for RHICE this may not be significant.) a is also expected to increase due to the effects of magnetic field in 
the plasma [lH|, further increasing the value of t. 

It is important to note that the initial magnetic field will enter the medium in the longitudinal direction as the 
medium is only about 1 fm thick in that direction at the initial stage. For this the relevant quantity is the penetration 
depth 8 ^ (fiauj)^^^^ where uj is the angular frequency of electromagnetic wave. Initial magnetic field, being a narrow 
pulse of time duration t ~ 0.2 fm (typically the width of Lorentz contracted Nuclei, for RHIC energies), can be taken 
to have w ~ 30 fm^^. This gives the penetration depth of order 3 fm. Thus, the picture of magnetic field diffusing 
through the entire region of the plasma with typical length scale of several fm, and lasting with high initial peak 
values for time scales of several fm, is self consistent. (Though, note that the time scale of the decay of the field being 
of same order as the system size makes the assumption of quasi-static field, and to that extent the assumption of ideal 
magnetohydrodynamics, only marginally valid.) As significant flow anisotropy develops also in the time scale of order 
3-4 fm [1^, it is then reasonable to assume that magnetic field can be taken to be almost constant for this duration 
of flow development. In conclusion, induced magnetic field, with similar magnitude as the peak value of the initial 
magnetic field pulse, decays slowly and, for time scales relevant for our model, can be taken to be approximately 
constant. 

With this discussion, we continue to investigate the effects of (an approximately constant) magnetic field on sound 
waves in QGP produced in RHICE. For the non-relativistic plasma the effect of magnetic field can be simply described 
in terms of three different waves (tI. [T7|. There is a fast magnetosonic wave which, for small magnetic fields, becomes 
the ordinary sound wave. Its velocity is given by = + v\ sir? 9, where va = Bo/y/Anp is the Alfven velocity 



and 9 is the angle between the wave vector and the magnetic field Bq. Cs is the sound velocity and p is the plasma 
density. Slow magnetosonic wave has velocity c?_ = v\ cos^ 9, and the Alfven wave velocity is given by c\ = w^os^ 9. 
These expressions for magnetosonic wave velocities are valid for the case when va « Cs- It was argued in [7| that 
fast magnetosonic waves lead to distortion of CMBR acoustic peaks while the slow magnetosonic waves may lead to 
long period modulation of the peaks. 

For the situation of plasma in RHICE one needs to consider the situation of relativistic magnetohydrodynamics. 
We will assume the applicability of magnetohydrodynamical description for the plasma at RHICE and restrict our 
consideration to wavelengths much larger than the Debye screening length. For the relativistic case, the expressions 
for the velocities of these waves are given below (we use natural units with c — 1). For the following discussion, we 
have followed ref. fl7{ . 



Alfven waves 

Phase velocity and group velocity Vg^ of these waves are 



Bo cos 6* Bo 

Vp/i = ^^n, Vgr = —=, (1) 
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where n = k/fc, k is the wave vector and 9 is the angle between n and the magnetic field Bq. — po^o + 
where h is the specific enthalpy (subscript denoting the background values) defined as [l7[ ph = pc^ + pe + P. Here, 



pe is the internal energy, pc is the rest mass energy and P is the pressure. For ultra-relativistic case we take ph — AP. 
Magnetosonic waves 

These are the waves which are relevant for our case of discussion of flow as they involve density perturbations. 
Phase velocities for these waves are given by 

= vph-^ = n(i[(po/io/wo)c? + v\]f/\\ -f <5cos2 9 ± af'^. (2) 
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FIG. 1: 

Typical situation expected in RHICE. Magnetic field points in y direction and the direction of the group velocity Vg 

obtained from n and t via Eq.(5). 



Here + and — signs correspond to the fast and slow magnetosonic waves respectively, 
speed, and S and a are defined below. 
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(cr here should not be confused with the conductivity discussed above.) For propagation of density perturbations, 
as relevant for the evolution of flow anisotropies, the relevant wave velocity is the group velocity for the magnetosonic 
waves, 
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[cr =F 2S{a ±{1 + S cos^ 0))] sin 9 cos t 
2(1 + (5 cos2 6* ±a)a 



(5) 



Here t — [(B0/-B0) x n] x n, and again the upper and lower signs (± or =f) correspond to the fast and the slow 
magnetosonic waves respectively. For a given magnetic field Bq, the direction of n can be varied to generate group 
velocities of these waves in different directions. Fig.l shows a typical situation of various vectors in Eq.(5) expected 
in RHICE. It is important to note that the direction of Vgr depends on the relative factors multiplying n and t in 
Eq.(5). This in turn depends on properties of the plasma like energy density. Thus due to the presence of spatial 
gradients in RHICE, even along a fixed azimuthal direction, we will expect the direction of Vg^ to keep varying with 
the radial distance. This can lead to the development of very complex flow patterns, possibly leading to generation 
of vorticity. 

For the strength of the magnetic field expected in RHICE, we use the estimates given in refs. d, [13] for Au-Au 
collisions at RHIC. For first estimates, we use a simple parametrization and take the magnetic field to be proportional 
to the impact parameter (with somewhat larger magnitudes than in ref. [T^] to illustrate the effects) . 
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(6) 



This corresponds to the maximum magnetic field of about 10^^ Tesla and the scale of 10 fm is taken from [l^ (where 
time integrated value of bBq is calculated, as coupling to charge particles is of interest there) as the approximate 
limiting value of the impact parameter b up to which proportionality with b is expected. The effect of magnetic 
field on wave propagation here comes from an effective magnetic pressure arising from the freezing of magnetic field 
lines in the plasma in the magnetohydrodynamic limit. Note that the importance of charges (quarks) here is only in 
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establishing the magnetohydrodynamic conditions such that the electric field in the comoving frame vanishes. There 
is no direct pressure generated by interaction of magnetic field with quarks. (Thus it is of no relevance that gluons, 
which do not interact with the magnetic field, dominate the energy density of the plasma.) The distortions of magnetic 
field lines in the presence of density perturbations cost energy leading to an extra contribution to pressure from the 
presence of magnetic field. This is what is responsible for increasing the effective sound speed as given above. 

We will assume a magnetic field with magnitude Bq given above which is uniform over the region of the plasma, 
and as we have discussed above, is approximately constant for the early stages of time scale of few fm. This early time 
duration is important for the evolution of flow anisotropics [l6| . The magnetic field has direction in the transverse 
plane, normal to the direction of the impact parameter. From the expression of Vg,. in Eq.(5) we see that an important 
factor is the ratio Bq /P where P is the pressure of the plasma which we take to be the quark-gluon plasma with two 
light flavors with the pressure given by P — ^n^T'^. Again, to illustrate the effects of the magnetic fleld on flow, we 
consider the situation at a lower value of the temperature T = 180 MeV. As the important effects occur for strong 
magnetic field which occurs for large impact parameter, a lower value of T (compared to what is expected in central 
collisions) may not be very unreasonable. We are considering the effects of magnetic field relevant for Au-Au collisions 
at RHIC on somewhat larger side of the estimates. However, for higher energy collisions, e. g. a t LHC, larger values 
of Bo (by an order of magnitude compared to RHIC) should be routinely expected, see e.g. flSj . 

We follow the procedure described in [l| for calculating w™* using HIJING [l^ . We start with the initial transverse 
energy density €tr distribution for Au-Au collision at 200 GeV/A center of mass energy from HIJING. For details, 
see ref. We assume that the hydrodynamic description becomes applicable by t = Tgg, which we take to be 1 fm 
and calculate the anisotropics in the fluctuations in the spatial extent i?((/)) at this stage, where i?(0) represents Ctr 
weighted average of the transverse radial coordinate in the angular bin at azimuthal coordinate (p. As emphasized 
above, angle (j) is taken in a lab fixed coordinate frame. We divide the region in 50 - 100 bins of azimuthal angle 0, 
and calculate the Fourier coefficients of the anisotropics in SR/R = {R— R{(j>))/R where R is the angular average of 
R{4'). Note that in this way we are representing all fluctuations essentially in terms of fluctuations in the boundary 
of the initial region. We use F„ to denote Fourier coefficients for these spatial anisotropics, and use w„ to denote 
nth Fourier coefficient of expected momentum anisotropy in Sp/p defined in the lab frame. We have generated events 
using HIJING and we present sample results for Au-Au collision at 200 GcV/A center of mass energy. In all the plots, 
the averages are taken over 1000 events. 

In I'l, the root mean square values v™^ of the fiow Fourier coefficients were obtained from spatial FnS simply by 
using proportionality factor of 0.2 (with a minus sign as Fn will be negative). We include the effect of magnetic field 
and the resulting angle dependent velocity of the (fast) magnetosonic wave by replacing this proportionality factor 
to 0.346 xvgr- Here the group velocity of the (fast) magnetosonic wave v^^ (Eq.(5)) changes with the angle in the 
event plane. The assumption here is that the flow coefficients are proportional to the sound velocity [iGj . (Ignoring 
that now the time scale of the development of flow may also vary with the azimuthal angle). The factor of 0.346 is 
chosen so that the proportionality constant becomes 0.2 for zero magnetic field case with the usual sound velocity 
cs = l/\/3. 

In Fig. 2 we show the plots of f™" for different values of impact parameter b with magnetic field given by Eq.(6) 
(Solid curves) and without magnetic field (dashed curves). These plots show that magnetic field can strongly affect 
values of v™^. 

We have also calculated average values of flow coefficients in the event plane. Fig. 3 shows the plots of the average 
flow coefficients w„ for 6 = 10 fm. Solid curve shows the plot with the presence of magnetic field and the dashed curve 
shows the plot in its absence. Note that both the curves approach almost zero value beyond n ~ 10 (as higher VnS 
lose correlation with the event plane) . 

It is important to study the effect of magnetic field on the elliptic flow V2 ■ In Fig. 4 we show the plot of the ratio 
of V2 with and without the magnetic field, i.e. V2{B)/v2{0) as a function of the impact parameter b. Elliptic flow 
itself varies with b and this ratio helps us in separating the effect of the magnetic field on V2 as the field varies with 
b. Note that this ratio becomes as large as 1.3 for 6 = 10 fm (for which the magnetic field takes its largest value 
in Eq.(6)). This strong enhancement in the value of V2 is important especially as it raises the interesting possibility 
whether a larger value of 77/s can be accommodated by RHIC data when proper accounts of magnetohydrodynamics 
are incorporated in flow calculations. 

Due to anisotropy of magnetosonic wave velocity, the radial flow itself will get modiflcd. To incorporate this effect 
we modulated the initial spatial proflle with a suitable weight factor proportional to Vg^- However, its effects were 
insignificant. It still remains a possibility that for much stronger magnetic fields non-trivial flow anisotropics may 
arise even with almost isotropic initial conditions (though, for central collisions one expects very small magnetic 
fields). Also, note that the direction of Vgr depends on the relative weights of factors multiplying n and t in Eq.(5) 
which depends on quantities like the plasma energy density etc. Thus due to the presence of spatial gradients of 
plasma density in RHICE, even along a given azimuthal direction in the transverse plane the direction of Vgr will keep 
changing. In particular, fluctuations of energy density will lead to fluctuations in Vg^ as well. Clearly due to such 
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FIG. 2: 

Plots of V™" for different values of impact parameters b. Solid (dashed) curves show the plots in the presence (absence) of 
magnetic field. (a),(b),(c),(d) correspond to the values of & = 10,8,6,2 fm respectively. 




FIG. 3: 

Plots of the average flow coefficients Vn for 6 = 10 fm. Solid curve shows the plot with the presence of magnetic field and the 

dashed curve shows the plot in its absence. 




FIG. 4: 

Plot of the ratio of V2 with and without the magnetic field, as a function of the impact parameter b (in fm). 
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phenomena one expects a complex pattern of flow, even possibly leading to vorticity, to develop in RHICE than just 
radial flow and flow anisotropics. We have also checked the effects of magnetic field on other features of the plots of 
^rms discussed in [![ , in particular on the acoustic peaks. The effects are similar to what is shown in Fig. 2 hence we 
do not show it here. It remains to be explored how the slow magnetosonic waves affect various features of the plots 
of v™^ in view of its proposed long period modulation of the acoustic peaks in CMBR. The role of Alfven waves also 
needs to be explored in the evolution of various fluctuations in RHICE (especially in view of their effects on CMBR 
[20j ) . Our approach in this work has been to present rough estimates of various possible effects of the presence of the 
magnetic field in RHICE. Detailed magnetohydrodynamical simulations are needed to probe these effects. Especially 
exciting will be the possibility of large effects at LHC energies where large magnetic fields (by almost an order of 
magnitude compared to RHIC) are expected [isj . 

We are very grateful to Abhishek Atreya, Anjishnu Sarkar, Uma Shankar Gupta, and Trilochan Bagarti for useful 
discussions. 
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